function PLoss = PLoss_Calculation( m_Node, Node_Num, m_G, m_B )
PLoss = 0.0;
for i = 1:Node_Num
    for j = 1:Node_Num
        f = m_Node(i).Angle - m_Node(j).Angle;
        PLoss = PLoss + m_Node(i).Vol * m_Node(j).Vol * ( m_G(i,j)*cos(f) + m_B(i,j)*sin(f) );
    end
end
